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We use CMB (WMAP and ACBAR), large scale structure (SDSS luminous red galaxies) and 
supernova (SNLS) data to constrain the possible contribution of CDM isocurvature modes to the 
primordial perturbation spectrum. We consider three different admixtures with adiabatic modes in a 
flat ACDM cosmology with no tensor modes: fixed correlations with a single spectral index; general 
correlations with a single spectral index; and general correlations with independent spectral indices 
for each mode. For fixed correlations, we verify the WMAP analysis for fully uncorrelated and anti- 
correlated modes, while for general correlations with a single index we find a small tightening of the 
constraint on the fractional contribution of isocurvature modes to the observed power over earlier 
work. For generally-correlated modes and independent spectral indices our results are quite different 
to previous work, needing a doubling of prior space for the isocurvature spectral index in order to 
explore adequately the region of high likelihood. Standard Markov-Chain Monte-Carlo techniques 
proved to be inadequate for this particular application; instead our results are obtained with nested 
sampling. We also use the Bayesian evidence, calculated simply in the nested-sampling algorithm, to 
compare models, finding the pure adiabatic model to be favoured over all our isocurvature models. 
This favouring is such that the logarithm of the Bayes Factor, ln_B < —2 for all models and lnfi < — 5 
in the cases of fully anti-correlated modes with a single spectral index (the curvaton scenario) and 
generally-correlated modes with a single spectral index. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

The most recent observations of the cosmic microwave 
background radiation (CMB) by the WMAP experiment 
[1] have been shown to be consistent with a now-standard 
cosmological model 2 , in which nearly scale-invariant, 
Gaussian, adiabatic primordial perturbations act as the 
seed of structure formation. This simple picture for the 
primordial perturbations is consistent with the expecta- 
tion of single-field inflation models. However, high en- 
ergy physics motivates several models for the early uni- 
verse that generate additional features in the primordial 
perturbation, such as non-Gaussianity and isocurvature 
modes, and constraining these features is an essential 
part of testing such models. In this paper, we focus on 
an isocurvature component. Although purely isocurva- 
ture primordial perturbations were ruled out by measure- 
ments determining the location of the first acoustic peak 
in the CMB angular power spectrum [3J, many theories 
predict a non-negligible contribution from isocurvature 
modes which cannot currently be ruled out by the data. 
An important example is models of inflation with multi- 
ple scalar fields, whereby fluctuations in field space that 
are perpendicular to the background trajectory describe 
isocurvature modes [H [5] . Well-studied special cases in- 
clude the curvaton scenario [5] and the axion model [TJ. 
Whether isocurvature fluctuations generated during in- 
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flation survive to the radiation-dominated epoch is de- 
pendent on how the fields decay at the end of inflation 
and any subsequent thermalisation. The best-motivated 
scenarios lead to the cold dark matter (CDM) isocurva- 
ture mode, in which the CDM is perturbed relative to 
the other species, and this is the only mode we consider 
in this work. 

Making use of some of the most recent cosmological 
data available - including measurements of the CMB 
power spectra by WMAP [8] and ACBAR [9] and the 
clustering of luminous red galaxies in the SDSS [10 , as 
well as somewhat older supernova data from the SNLS 
[TTj - we readdress the question of constraints on CDM 
isocurvature contributions. This partially revisits a pre- 
vious analysis of the WMAP three-year data by Bean 
et al. [12] who suggest that improvements in large-scale 
measurements of polarization in WMAP five-year data 
should improve the constraints on isocurvature modes 
(based on the analysis of [IS])- In considering the sim- 
plest models (models containing fixed correlations of adi- 
abatic and CDM isocurvature modes) we do indeed see 
an improvement, and in doing so verify the results of the 
WMAP team's own analysis [Hj. We also see a slight 
improvement for models with generally-correlated adia- 
batic and CDM isocurvature modes with a single spectral 
index. However, our result for models with generally- 
correlated adiabatic and CDM isocurvature modes, with 
independently varying spectral indices for each mode, are 
very different from the published results in [12] - a dif- 
ference which presented itself consistently when we em- 
ployed the same (earlier) data set as in their analysis. In 
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correspondence with the authors of [T5] , we can attribute 
the difference to a bug in their code used to compute the 
cross power spectrum in models with correlated adiabatic 
and isocurvature modes and unequal spectral indices. 

In our analysis, we experienced some difficulties in 
applying standard Markov-Chain Monte-Carlo (MCMC) 
techniques to explore the parameter space of our most 
general models, particularly in relation to proposal dis- 
tributions and convergence statistics (see Sec. VBl. 



To overcome these problems we use a nested-sampling 
method (MultiNest) [TSl US], which we use to obtain 
results throughout this paper except where specifically 
stated otherwise. The agreement with the results of the 
WMAP5 analysis [14 , as well as consistency with some 
of the results of Bean et al. [T3] we believe provides 
appropriate validation of the nested-sampling method in 
this application. There is also some degree of consistency 
to be found with the earlier work of Beltran et al. [17] . 
our results showing a similar favouring of models with 
a high value for the isocurvature spectral index (> 3) 
when it is allowed to vary independently of the adiabatic 
spectral index. Keskitalo et al. [TS] also find a sim- 
ilarly large isocurvature spectral index to be favoured. 
These particular comparisons are tenuous given the use 
of slightly different paramcterizations to our own, but it 
is undoubtedly interesting to note these preferred mod- 
els show a similar excess of power on small scales in the 
matter power spectrum. 

A further benefit of the nested-sampling technique is 
that we can easily compute the Bayesian evidence, allow- 
ing us to perform a simple model selection analysis. By 
considering the Bayes' factor (evidence ratio) of two mod- 
els, we may use the Jeffreys' scale [19 to say which model 
is favoured by the data, and to what extent this favour- 
ing is significant. We are thus able to compare models 
including isocurvature modes to each other and to a pure 
adiabatic model. Trotta [20] provides a full review of in- 
terpreting the evidence and in this work we use the weak, 
moderate and strong definitions of the Jeffreys scale for 
the logarithm of the Bayes factor given therein. Some 
previous work in isocurvature model selection has been 
inconclusive [2T] . Other work [23] formulated the analysis 
such that a prior was placed directly on the isocurvature 
fraction. In the latter case, model comparison suggested 
a strong disfavouring of models that allowed isocurvature 
fractions larger than unity. 

In this work we have also carried out, for the first time 
in such complex isocurvature models, an analysis that 
includes the effects of weak gravitational lensing of the 
CMB by large scale structure (see [23] for a recent re- 
view) . Note that this required some modifications to the 
implementation of lensing in the CAMB code that we 
employ [55], since the total CMB angular power spectra 
must be lensed by the total lensing-potential power spec- 
trum [35] (i.e. including the correct isocurvature contri- 
bution). However, the effect of lensing is minimal: in our 
model of generally-correlated modes with independently 
varying spectral indices we see essentially no differences 



in the posterior distributions with and without lensing 
included. This suggests that, for these models, including 
lensing of the CMB in the analysis is unnecessary with 
current data. To avoid complicating our presentation at 
no net effect, we only present results obtained without 
lensing here. 

The structure of this paper is as follows. In Sec. [IT] 
we introduce the parameterizations used in our analysis 
and our background cosmology. In Sec. |III| we introduce 
the data sets, priors and sampling methodology we use 
to generate constraints. In Sec. |IV| we present our re- 
sults for models with an admixture of adiabatic and CDM 
isocurvature modes and fixed correlations between these 
modes. We extend this to general correlations in Sec. \V\ 
and discuss our findings and their possible implications 
in Seclyll 



II. PARAMETERIZATIONS 

In the standard cosmological model the primordial per- 
turbation is purely adiabatic: all species are perturbed 
similarly so that the gauge-invariant entropy perturba- 
tion 
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vanishes for all i and j. The regular mode of such pertur- 
bations is characterized by a single field which we take to 
be the comoving-gauge curvature perturbation 1Z. The 
adiabatic condition is preserved on super-Hubble scales 
and 1Z is conserved. In contrast, isocurvature perturba- 
tions have compensating density perturbations that pro- 
duce no net curvature at early times, 1Z — 0. In this pa- 
per we consider only the CDM isocurvature mode which 
is characterized by the non-zero S = S cr between the 
CDM and the radiation. Our sign conventions are such 
that, on large angular scales, the Sachs- Wolfe contribu- 
tion to the temperature anisotropies is |25j 
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Note that this convention for 1Z is opposite to that some- 
times employed, e.g. [2"5] . 

The most general model we consider contains an ad- 
mixture of the adiabatic and CDM isocurvature mode 
with general correlations. We characterize the relevant 
two-point correlation functions by power-laws: 
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where A are amplitude parameters, n are spectral in- 
dices and fco is a pivot value (taken throughout this work 
to be 0.05 Mpc -1 ). Subscripts A and / label adiabatic 
and isocurvature modes respectively and we take hai — 
(jia + nj)/2. This choice for hai makes the adiabatic- 
isocurvature correlation coefficient independent of scale, 
and so positive-definiteness of the matrix Aij of ampli- 
tude parameters is sufficient to ensure the same for the 
primordial spectra at all scales. Note, however, that in 
general two-field inflation models, there is no such sim- 
ple relation between the three spectral indices [27], but 
allowing uai to vary independently was shown to have 
little effect in previous analyses [17 . 

For photon transfer functions Op^^fc), where i = A 
or / labels the mode and X denotes the CMB observable 
(T or E here), we define a c\ XY ^ as 
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We may then write the total CMB power spectra, Ci, as 

C t = A AA Cf A + AuC 1 / + 2A AI Cf , (5) 

for all spectra XY. This is commonly parametrized by 
a and (3 such that 
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where A — Aaa + Ah is an overall amplitude parameter. 
In this parameterization a, lying in the range [0,1], quan- 
tifies the amplitude of isocurvature contribution with 
a/(l — a) = Au/Aaa- The parameter 
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with a range [—1, 1], is the correlation coefficient between 
the two modes. We will refer to this parameterization as 
the a/3-parameterization. It should be noted that the 
definition of (3 is consistent with that given in [2] . Note 
also that Aij, C\\ A and a all depend on the pivot scale 
ko, but (3 is independent of fco- We also construct the 
total matter power spectrum from the primordial spec- 
tra using the matter transfer functions for adiabatic and 
isocurvature modes from CAMB. 

An alternative parameterization, motivated by the de- 
sire to include more than two perturbation modes, is 
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where the C^ Y ^ have been normalised by the 
temperature-anisotropy power Pi = ^1=2° l (2^+l)cj TT ' ) " 
such that C 1 / = C 1 / / yJPjPj. In this case, ^JPiPjAij cx 
Zij with Zij a positive-definite matrix normalised such 
that 
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= 1. We sample the Zij uniformly 
from the surface of the unit 2-sphere in the basis 
(zaa, zn, V^zai), rejecting points for which the result- 
ing Zij is not positive definite. Rather than sampling 
in B, we sample in the total (TT) CMB power (up to 

,1000 



1000), B = Ee=2 ( 2i + 1)PP so that 
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We will refer to this parameterization as the z- 
parameterization. It has the virtue of being more closely 
linked to the observable isocurvature contribution than 
a-/3; in particular, z^ and B are independent of the pivot 
scale k n . Furthermore, the correlation coefficient f3 is 
unconstrained by the data near a = and a = 1 [see 
Eq. ([6|] but the z-parameterization does not suffer from 
this problem. In the z-parameterization, the amount of 
isocurvature can be quantified by the derived parame- 
ter r iso = Zi S o/(2:i SO + z AA ) where z iso = y/l - z\ A [T2]. 
This represents the fractional contribution of isocurva- 
ture modes to the observed TT power at I < 1000. 

We shall use both parameterizations in this work, to 
enable comparison with previous work and to explore the 
dependence of isocurvature constraints on the choice of 
prior. In our case of a single isocurvature mode mixed 
with the adiabatic mode the relations between (a, [3) and 
the Zij are 
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In this work we use a background cosmology that is a 
flat ACDM model. We parametrize the model by the den- 
sity parameters, fi&/i 2 , £l c h 2 , Qa, and the optical depth r 
in addition to A or B, two isocurvature contribution pa- 
rameters [either (a, (3) or the amplitude parameters for 
mapping to z^] and spectral indices tia and rij. In using 
WMAP and ACBAR data we marginalise over the am- 
plitude of a template for the thermal Sunyaev-Zel'dovich 
(SZ) effect since this is non-negligible at the scales now 
covered by these experiments. We ignore any dependence 
of the SZ template on the isocurvature contribution. This 
should be harmless for models with nj « 1, since the con- 
tribution of isocurvature modes to the late-time small- 
scale power is then small, but may be worth revisiting 
in future analyses of models with very blue isocurvature 
spectra when better quality data is available on small 
scales. 
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TABLE I: The basic flat prior ranges for our sampling pa- 
rameters. Priors on flth 2 ,Q c h 2 and Q.a are affected by a BBN 
prior and Ho prior giving rise to the effective priors plotted 
m Fig. [I] Where we constrain ni = ua we use the ua prior. 



The difference in the priors of the z-parameterization 
and the a/3-parameterization can be most clearly seen 
in Fig. [2j in which we have used Eqs. (Ill and (12) to 



Parameter 



Range 
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r 
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[0.005,0.1] 
[0.01,0.99] 

[0,0.9] 
[0.01,0.8] 
[0.8, 1.2] 

[0,6] 
[2.6,4.2] 
[0,1] 



III. DATA SETS AND PRIORS 

We find constraints using CMB (WMAP5 [5] and 
ACBAR [9J, large scale structure (SDSS luminous red 
galaxies [10]) and supernova (SNLS [11]) data. We place 
a cut-off on the SDSS data, discarding points above 
k = 0.08/iMpc -1 to exclude scales where non-linear clus- 
tering is significant. We also impose a cut-off in the 
ACBAR data, ignoring £ > 2000 to reduce the impact on 
the analysis of foregrounds, such as dusty galaxies, and 
secondary anisotropics such as the SZ effect. A Gaussian 
prior of il^h 2 = 0.022 ± 0.006 is included to account for 
big bang nucleosynthesis (BBN) estimates of the baryon- 
to-photon ratio. We include no age prior but do require 
Hq to lie between 40 and 100 kms -1 Mpc -1 . Subject to 
these constraints, and those outlined in Sec.[ll] the priors 
we use for sampling are shown in Table [I] 

Figure [TJ shows the set of Id marginalised priors for the 
z-parameterization, omitting parameters with flat priors, 
obtained from MCMC runs of CosmoMC [29] with no 
data. We also include in this figure the effective prior 
on rj so . The requirements that Zij is a positive-definite 
matrix and ^ zfj = 1 are clearly visible in their Id 
marginals. The effect of our BBN and Hq priors can 
also be seen in the distributions for Q(,ft. 2 , £l c h 2 and 
The priors on these parameters are the same in the a(3- 
parameterization, which has flat priors on a and (3 over 
the ranges [0, 1] and [—1,1] respectively. The effective 
prior on r iso can be calculated in parametric form, 

Pr(nso) oc sin 0(1 + sin 20) yn — 2 sin - 1 (tan 0/2)} , (13) 

where r- lso — sin 0/ (sin + cos0); this is also plotted in 
Fig. [T] The prior favours roughly equal contributions 
from adiabatic and isocurvature modes to the observed 
power, and goes to zero at r; so = and r[ so — 1. Note 
that if the data favours values around rj so = 0, the van- 
ishing of the prior at rj so = will force a peak in the pos- 
terior away from zero. This deficiency in the parameter 
riso = could easily be remedied by defining alternative 
measures of the isocurvature fraction, although we do not 
pursue this here |36] . Related prior issues have also been 
noted in Ref. |2"51. 



map the priors on to a and (i. In the top row we 
show how flat priors on a and (3 translate to a prior 
distribution in the coefficients of C 1 / and Cf 1 [a and 
2/3y / a(l — a)}. In this plane, the prior volume is en- 
closed in the ellipse [/3y/a(l - a)} 2 + (a - 1/2) 2 = 1/4 
and the prior density is oc 1/ \J a(l — a). In the bottom 
row of Fig. [2] we see how the effective Id marginalised 
priors on a and (3 in the z-parameterization are very 
non-uniform, with a = 1 (i.e. pure isocurvature) strongly 
favoured. This preference follows from Pa 2> Pi for mod- 
els with tia ~ nj since, to generate equal observed power 
as favoured by the z-parameterization, the primordial 
isocurvature power must be much larger than the adi- 
abatic. We find that the strong favouring for a ~ 1 in 
the z-parameterization prior is overwhelmed by the data 
in all cases. 

In a complete analysis, we would explore the posterior 
for each of the four data sets separately, before carrying 
out the exploration with the full data set. This would 
then allow us to check the consistency within each model 
of the data sets, and would also allow us to assess the con- 
straining power of each data set upon the model via the 
Kullback-Leibler divergence Z?kl [H3I3I]- Unfortunately, 
however, it is currently computationally prohibitive to 
repeat our analysis four times over, given the long time 
that it takes to calculate the spectra involved. Previous 
analyses of simple models including isocurvature [TH Q3] 
indicate consistency between WMAP and data sets in- 
cluding WMAP, supernova and matter power spectrum 
data. Additionally, Beltran et al. [T7] suggest that mat- 
ter power spectrum constraints on complex isocurvature 
models will only be significant for accurate measurements 
on small scales, which we do not have with our data set. 
We believe it is reasonable to assume therefore that the 
main constraining power on our more complex models 
comes from WMAP with other data improving our con- 
straints in a consistent manner. 

The posteriors are explored using the MultiNest nested 
sampler |16] incorporated into CosmoMC [29 , which em- 
ploys CAMB [21] to generate Cg and matter power spec- 
tra. We use 200 live points and a sampling efficiency 
of 1.0 which we believe is suitable for giving an indi- 
cation of which model the data favours in this appli- 
cation. The nested sampler was found to be necessary 
for generally-correlated models with independent spec- 
tral indices since, in these models, traditional MCMC 
techniques tended to get stuck in certain areas of the 
likelihood distribution, resulting in either the oversam- 
pling of areas where the likelihood was relatively low or 
the non-discovery of other areas of higher likelihood. No- 
tably, the Gelman-Rubin convergence statistics in cither 
case suggested erroneously that the chains were reason- 
ably converged [37]. This was likely due to the need to 
make the proposal distributions too narrow in order to 
explore the distributions with reasonable computational 
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FIG. 1: The non-flat Id marginalised prior distributions on the sample parameters and the effective prior on the derived 
parameter ri so in the ^-parameterization. These priors were obtained from MCMC runs of CosmoMC with no data. The red 
line in the r- lso plot is the analytic form of the prior, Eq. ( |13[ ), showing clearly how the prior goes to zero for r- lso = and 1. 
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FIG. 2: One-dimensional marginalised priors on a and j3 and 2d prior surfaces in a and 2/?ya(l — a) for the two parameter- 
izations. In the top row we show how the flat priors on a and j3 in the a/3-parameterization map to a prior distribution in the 
coefficients of Cl 1 and Cf 1 . The combination of sampling priors in the z-parameterization is mapped to effective priors on a 



and (3 by Eqs. (Ill and ( 12 1 in the bottom row, with the resulting prior distribution in a and 2/3 1 /a(l - a) displayed at the 
end of the row. The prior distributions were generated from MCMC runs of CosmoMC with no data. 



efficiency. The nested sampler, which converges on the 
posterior distribution from an initial uniform sampling 
of the prior, does not have these problems, and is used 
to obtain results for all models throughout this paper 
for consistency. Details of the difficulties we encountered 
with MCMC methods are discussed further in Sec. IV Bl 



IV. UNCORRELATED AND 
ANTI-CORRELATED MODES WITH A SINGLE 
SPECTRAL INDEX 

We first consider uncorrelated and anti-correlated 
modes in the a/3-parametcrization ((3 = and (3 = — 1 re- 
spectively) with the added condition that rij = = n. 
These two models could be produced by the axion [7] 
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FIG. 3: The Id marginalised posterior distributions for a in 
the cases j3 = and /3 = — 1 and 2d posterior contours at 
68% and 95% confidence level for a and the spectral index n. 
The 95% confidence limits on a are a < 0.15 for /3 = and 
a < 0.011 for /3 = -1. 



and curvaton [6] scenarios respectively, and were consid- 
ered in the light of the WMAP5 release in [14] (though 
with a pivot ko = 0.002 Mpc -1 and nj — 1). Here, for 
direct comparison with |14j we use only the most recent 
WMAP TT, TE and EE spectra to generate posteriors. 
It follows from Eq. Q that adding isocurvature modes 
boosts the CMB power on large scales for both (3 = 
and (3 = — 1. 

In the case (3 = 0, we find confidence limits at the 95% 
level to be a < 0.15. For (3 — — 1 we find the 95% limit 
to be a < 0.011. As expected, the constraint on a is 
tighter for anti-correlated isocurvature modes since they 
have a proportionately larger effect on the large- angle 
CMB spectra. Both our results are in very good agree- 
ment with [2] which finds 0.16 and 0.011 respectively. 
Furthermore, the shapes of the marginalised posteriors 
for the two models shown in Fig. [3] are very similar to 
their counterparts in |14j with a clear degeneracy visible 
in the contours between a and n. With ha = n.j, the 
observed effect of isocurvature modes is limited to large 
scales and can be offset by increasing n. The inclusion 
of our other data sets would help to break this degener- 
acy, either by directly improving the measurement of n 
or indirectly by improving the measurement of the mat- 
ter density with which ua is partly degenerate in current 
CMB data. An example of the latter is provided in [TJ] 
where constraints on a tightened considerably once dis- 
tance information (supernovae and baryon acoustic oscil- 
lation data) was added. A full list of confidence limits 
may be found in Table [TT] Most parameters are stable 
to the inclusion of isocurvature modes, with mean val- 



TABLE II: Means and 68% confidence limits or 95% one-tail 
confidence limits for models with fixed correlations of adia- 
batic and isocurvature modes. 



Parameter Pure Adiabatic 



Q b h 2 
Q c h 2 



in 



n 
[10 10 



u.uzzj_ 0006 
0.109 ±0.006 
0.744 ±0.029 
0.0887lHi 5 5 t 
0.960 ±0.013 
3.06 ± 0.04 



(3 = 



(3 = -\ 



0.0231 ±0.0008 
0.105 ±0.007 
0.768 ± 0.032 
0-0880l°;gi» 
0.980t° °! 

< 0.152 



0.0229 ± 0.0006 
0.103 ±0.006 
0.779 ±0.027 

0.0879tS:SiE 
0.980 ±0.018 
3.04 ±0.04 
< 0.108 



ues shifting by less than la. The largest effect is seen in 
n, where including isocurvature modes shifts the mean 
upwards by 1.5c, due to the degeneracy with a. We 
regard our WMAP-only analysis as adequately demon- 
strating the robustness of the nested-sampling method 
in the current application. 

The logarithm of the Bayes' factor B between the axion 
model and the pure adiabatic model is found to be lni? = 
—2.13 ±0.43, favouring the pure adiabatic model moder- 
ately on the Jeffreys' scale. Between the curvaton model 
and pure adiabatic model, we find ln_B = —5.97 ± 0.45, 
which can be interpreted on the Jeffreys' scale as strong 
evidence for the pure adiabatic model over the curvaton, 
based on WMAP data alone. 



V. GENERALLY-CORRELATED MODES 

More general correlations between the adiabatic and 
isocurvature modes arise naturally in two-field inflation 
when the background trajectory is curved in field space 
after horizon crossing (e.g. |27]V We carry out our anal- 
ysis of such generally-correlated modes only in the z- 
parameterization, with a and /?, where they appear, hav- 



ing been mapped from the z-matrix using Eqs. ( 11 ) and 
(12 1. We consider two cases, ua — nj and ua and nj 



varying independently, and use all four datasets in both 
cases. 



A. Single spectral index 

The (marginalised) 68% and 95% contours in the zai- 
zn plane are plotted in Fig. [4] for the case ha = nj. 
These are very similar in shape to the results in |12] 
(Fig. 2 there) which were obtained with earlier data. 
Note that we have applied the prior constraint of positive- 
definiteness, z AI + (zjj — 1 /2) 2 < 1 /4, directly in the plot 
which gives the lower hard edge to the contours. We find 
a slight tightening in the marginalised confidence lim- 
its over those in [12], as expected from the more recent 
data: zn < 0.08 and zaa > 0.992 (at 95% confidence) 
compared to z H < 0.09 and z A a > 0.989 from [TJ . The 
mean and 68% confidence limits zai = 0.03 ± 0.03 also 
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FIG. 4: Contours at the 68% and 95% confidence level in the 
zai-Zii plane (left) and the Id marginalised posterior distri- 
bution (black) and prior distribution (red) for ri so (right) in 
the case of generally-correlated modes with a single spectral 
index. 




FIG. 5: Contours at the 68% and 95% confidence level show- 
ing the effect of independently varying the isocurvature spec- 
tral index ni on the correlation coefficient f3 (left) and the adi- 
abatic spectral index n A (right). Note that we have sampled 
in the ^-parameterization and mapped to ft using Eq. (|12[). 



TABLE III: Means and 65% confidence limits or 95% one- 
tail confidence limits for models with general correlations of 
adiabatic and isocurvature modes. 



Parameter Pure Adiabatic 



ni = n A 



Qch 2 
Oa 

T 
UA 

m 

ZAA 
ZAI 
Zll 



n 0226 +IMmu(i 
u.uzzo_ 0005 

0.108 ±0.004 

0.755 ±0.020 

0.0884^.^8 
0.964 ±0.013 



0.155 ±0.002 



0.0231 ±0.0007 
0.106 ±0.004 
0.761 ±0.022 
0.08681™ 
0.976 ± 0.020 

0.155 ±0.002 

> 0.992 
n fi9c;9+ 0303 

u.uzoz_ 02g8 

< 0.0783 

< 0.113 



n-i ^ n A 
u.uzoo_ 0010 
0.108±g;gg| 
0.766 ±0.021 

+0.0173 
-0.0164 



0.0869 
0.956 ±0.015 

q yo-t-0.38 
J - ' "-0.08 

0.155 ±0.002 

> 0.998 
0.00171S 
0.0156™ 

3+0.0044 
-0.0105 



0.03097 



improve on the previous values 0.06 ± 0.07. A full list 
of confidence limits can be found in Table IIIIl With the 
inclusion of all four datasets, all parameters are stable 
to the inclusion of isocurvature modes with mean values 
shifting by less than la. 

In right-hand panel of Fig. [4] we plot the Id 
marginalised distribution for r lso . The 95% upper con- 
fidence limit for this distribution is r- lso < 0.11 (0.13 in 
|12j). Note that the data clearly overwhelms the prior in 
limiting r- lso (to the extent that the first bin of the his- 
togram contains enough data to give the impression that 
the posterior distribution is non-zero for ri SO = 0) and 
also that a peak has appeared near ri SO = 0. However, as 
the distribution has only a one-tail confidence limit it is 
not a concern to what extent this peak is an artefact of 
the prior. 

The logarithm of the Bayes factor between the 
generally-correlated isocurvature model with tia = ni 
and the pure adiabatic model is h\B = — 6.84± 0.48. At 
this magnitude the Jeffreys' scale suggests the evidence 
is very strong in favouring the adiabatic model over the 
isocurvature model. 



B. Independent spectral indices 

When we allow the spectral indices to vary indepen- 
dently, our results are quite different to the earlier, erro- 
neous, analysis in [T2]. (As noted in Sec. |IJ the results 
of Sec. IV of [12] are affected by a bug in the code the 
authors used to generate Cf 1 for tia ^ nj.) For exam- 
ple, in Fig. [5] we show the 68% and 95% contours in the 
nj~P and ni-riA planes which can be compared directly 
to Fig. 3 in [12]. The most notable difference is that we 
found it necessary to double the extent of the n/ prior 
(to [0, 6]) to explore adequately the regions of high like- 
lihood. We obtain n/-/3 contours that are largely sym- 
metric about (3 — (no correlation) and so there is no 
preference for the anti-correlation reported in |12j . More- 
over, and significantly, we find that the detection of a red 
adiabatic spectral index tia < 1 is stable to the addition 
of CDM isocurvature modes (for models with no grav- 
itational waves). We find the 95% one-tail confidence 
limit zaa > 0.998 (0.996 in [T5]) and the mean value 
and 68% confidence limits zai — 0.002 ± 0.023 (com- 
pared to —0.02 ± 0.01). The 95% upper confidence limit 
on zu < 0.03 (0.05 in [T3]) is now accompanied by a 
lower limit of 0.005. Additional confidence limits may 
be found in Table |III| For the parameters describing the 
background cosmology, only the baryon density shifts sig- 
nificantly (upwards by 1.8<r for the mean) from the pure 
adiabatic model. 

It is in the particular case of independently varying 
spectral indices that the advantages of using the nested 
sampler become apparent. If we use the now-standard 
MCMC method (CosmoMC g^) in this application, it 
is computationally prohibitive to explore the relatively 
large volume of prior space with low likelihood with an 
(initially) broad proposal distribution. By reducing the 
width of the proposal distribution we are able to pass 
the initial "burn-in" phase, but this results in the chains 
becoming "stuck" , oversampling regions of the posterior 
and thus potentially distorting one's conclusions. Using 
an initial run to compute a covariance matrix for subse- 
quent runs also exhibits the same problem. 
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FIG. 6: One-dimensional marginalised posterior distributions 
for n/ (top) and a (bottom). The left-hand panels are ob- 
tained with the nested-sampling method; the right-hand pan- 
els show a set of MCMC chains which are converged (the 
Gelman-Rubin convergence statistic R = 0.06 for both m 
and a). 
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FIG. 7: One-dimensional marginalised posterior distributions 
for ni for each individual MCMC chain in the set used in the 
right-hand panels of Fig. [6] Though showing clear discrepan- 
cies, these chains may be considered to be converged by their 
Gelman-Rubin convergence statistic, R = 0.06. 



The phenomenon is particularly visible in two pa- 
rameters, 71/ and the derived a, and is illustrated in 
Fig. [6] The posterior distributions for these parameters 
obtained with nested sampling show a small rise around 
a = = nj , though with a small integrated weight com- 
pared to the main peak. The set of MCMC chains plotted 
in Fig. [6j appears to have become rather stuck around 
a = = ri/, oversampling those regions and distort- 
ing the likelihood inferred from the chains even though 
the chain has supposedly converged - the ratio R of the 
variance of the chain mean to the mean of the chain vari- 
ance has R < 0.07 for all sampling parameters [35]. In 
Fig. [7] the marginalised nj distributions for each of the 
four chains in the set are plotted separately, showing clear 
disparity between two chains that have spent much time 
near 71/ = and two chains that spent more time in the 
main peak. 

In Fig. [8] we plot the 68% and 95% contours in the Zai- 
zn plane, and the marginalised Id posterior for r; so in the 
case of independent variations in jia and ti/. Compared 
to Fig. [4] which enforced ua = 71/ , zjj is much more 
tightly constrained. In addition, the putative peak in the 
marginal distribution for r; so with Ua — ni is sharpened 
considerably around ri SO = 0.02. The 95% confidence 
limits are 0.01 < n so < 0.06 (compared to the one-tail 
upper limit of 0.10 in [H]). The extent to which this peak 
is driven by the prior is difficult to quantify. However, 
some insight may be gained by ignoring the first bin used 
to compute the histogram for ri SO (containing the prob- 
lematic zero in the prior) and calculating the ratio of the 
peak likelihood to the likelihood at the second bin. We 



find that this ratio is ~ 8 times larger than the ratio of 
the prior values at the same points. This would suggest 
that the peak is not driven by the prior and appears to 
be a real feature of the data. 

The best-fit model has a very blue isocurvature spec- 
trum, 71/ = 4.19 (hence the need for a wide prior range). 
The best-fit model improves the log-likelihood over the 
pure adiabatic model by five units with three additional 
degrees of freedom. However, the logarithm of the Bayes 
factor between this model and the pure adiabatic model 
is InB = —3.82 ± 0.49, which would normally be con- 
sidered to be moderate to strong evidence for the pure 
adiabatic model over this isocurvature model. The im- 
provement in fit is overwhelmed by the large increase in 
prior volume in the comparison of the evidence. When 
considering the evidence of the isocurvature model, how- 
ever, we should be cautious as we have little physical 
justification for our wide prior on 71/ . Given that the 
evidence may be thought of as the average of the like- 
lihood over the prior, we can see from the 71/ posterior 
(Fig. [6]) that if our prior for this parameter was instead 
flat over the narrower range [3, 5] the evidence would im- 
prove, although we calculate this improvement to be just 
0.97 units of log evidence and therefore the model is still 
moderately disfavoured against the pure adiabatic model. 
If we take the prior on 71/ to be flat over the range [0, 3] 
(as in [IS]), the evidence value becomes much lower as 
we would expect, worsening by 1.48 units of log evidence 
and suggesting the pure adiabatic model to be strongly 
favoured by the data in this case. 

There are some parallels to be drawn between our 
results and the results of Beltran et al. [T7] who 
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FIG. 8: 68% and 95% contours in the Zai-Zii plane (left) and 
the marginalised posterior distribution (black) for ri so (right) 
in the case that ua and ni are allowed to vary independently. 
The effective Id marginalised prior on r[ so is included (red) 
in the right-hand plot. 



used a similar model to what we have called the a(3- 
parameterization. In the case of independently varying 
spectral indices, they found the likelihood peaked for high 
(> 2) isocurvature spectral tilt. These favoured mod- 
els exhibited an excess of matter fluctuations on small 
scales. Though we use the z-parameterization, we have 
found nj = 3.78 ± 0.76 where the errors are la. In 
Fig. [9] we plot the CMB and matter power spectra for 
the best-fit model, which has nj — 4.19 and r lso — 0.019, 
along with the same for the best-fit pure adiabatic model. 
These models may be distinguished with future CMB 
data at high £ 7 though we note that this will require an 
excellent understanding of the instrumental beams and, 
for the temperature, careful accounting for secondary 
anisotropics and point-source contamination. (The am- 
plitude of the thermal SZ template in the best-fit general 
model is less than one half that in the best-fit pure adia- 
batic model. It may be worth investigating the potential 
correlation between this parameter and the isocurvature 
spectral tilt when better data becomes available on small 
scales.) The differences in the matter power spectrum 
are more significant, with the model including isocurva- 
ture modes having a clear excess of power on small scales 
over the pure adiabatic model. Later work by Beltran et 
al. [35] used the Lyman-a forest to constrain further the 
matter power spectrum at small scales; they found this 
squeezed the likelihood in nj to a peak around m = 2. 
Future Lyman-a analyses may have the same effect on 
the constraints presented here. 



VI. CONCLUSION 

The contribution of a CDM isocurvature mode to 
the spectrum of primordial perturbations has been con- 
strained in flat, ACDM models with no tensor perturba- 
tions using data from the CMB (WMAP and ACBAR), 
large scale structure (SDSS luminous red galaxies), and 
supernova (SNLS). Using the nested-sampling technique, 
we have obtained parameter constraints and Bayesian ev- 
idence values for model comparison with the pure adia- 
batic model. 
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FIG. 9: Spectra of our best-fit model of generally-correlated 
isocurvature modes with ni = 4.19 (black) and the best-fit 
pure adiabatic model (red dashed): TT, TE and EE spec- 
tra (top); and matter power spectrum with SDSS luminous 
red galaxies data (bottom; the black dashed line indicates 
where we placed a cut-off on the SDSS data). The model 
with isocurvature has ri so = 0.019 and exhibits a clear excess 
of matter clustering on scales k > 0.1/iMpc -1 in the matter 
power spectrum. 



We have considered models with fixed (i.e. axion or 
curvaton-like) and generally-correlated modes (e.g. from 
multi-field inflation). We find all models to be dis- 
favoured when compared to the pure adiabatic model 
using the Jeffreys scale, and with all but the axion- 
like model (uncorrelated isocurvature perturbations with 
ua — nj) strongly disfavoured on the Jeffreys' scale. Our 
constraints for models with fixed correlations of adiabatic 
and isocurvature modes are consistent with the recent 
WMAP analysis [14] , while there is a small improvement 
in the constraints on generally-correlated models with a 
single spectral index over the earlier analysis in |12j due 
to improvements in the data. In the latter case we find 
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the fractional isocurvature contribution to the observed 
power, ri so < 0.11 at 95% confidence. Our results for 
generally-correlated modes with independently varying 
spectral indices correct those in [TJ]. In these models, we 
find a small island of probability at very low isocurva- 
ture fraction, and a main peak with a small contribution 
from a very blue isocurvature spectrum (rij = 4.19 and 
fiso = 0.019 in the best-fit model). Navigating this com- 
plex space proved problematic for MCMC techniques but 
nested sampling avoids such problems. With nested sam- 
pling, we find the fractional (observed) isocurvature con- 
tribution 0.01 < rise < 0.06 at 95% confidence, but cau- 
tion that interpreting the lower limit is not straightfor- 
ward since the prior density vanishes at r iso = 0. We also 
find that a red adiabatic auto-spectrum, tia < 1, is ro- 
bust to the addition of general CDM isocurvature modes. 
The best-fitting model has an improved log-likelihood 
by six units compared to the adiabatic model but the 
larger prior volume negates this in the evidence which 
strongly favours the adiabatic model. We caution that 
the isocurvature spectral tilt is poorly constrained by the 
data, compared to the adiabatic tilt, and our (large) prior 
range is not strongly physically-motivated. However, the 
strong preference for adiabaticity appears stable to sim- 
ple changes in the prior. 

In our most general model we find there to be a strong 



preference for a very blue isocurvature spectral tilt. Such 
tilts exhibit a large excess on small scales in the matter 
power spectrum over the pure adiabatic model. While 
future constraints on the isocurvature contribution will 
undoubtedly be improved by better data from Planck |33j 
and other CMB experiments, we expect that the small- 
scale matter power spectrum will also play a crucial role. 



Acknowledgments 

This work was performed largely using the Dar- 
win Supercomputer of the University of Cam- 
bridge High Performance Computing Service 
( http://www.hpc.cam. ac.uk/] ), provided by Dell Inc. 
using Strategic Research Infrastructure Funding from 
the Higher Education Funding Council for England. We 
thank Michael Bridges, Farhan Feroz and Andy Taylor 
for helpful discussions. We are particularly grateful 
to Jo Dunkley for correspondence over the differences 
between our analysis of generally-correlated models 
and that in [15], and also for her help with alternative 
convergence tests for the MCMC chains. IS is supported 
by an STFC studentship. 



[1] M. R. Nolta et al. (WMAP) (2008), 0803.0593. 

[2] J. Dunkley et al. (WMAP) (2008), 0803.0586. 

[3] K. Enqvist, H. Kurki-Suonio, and J. Valiviita, Phys. Rev. 

D62, 103003 (2000), astro-ph/0006429. 
[4] D. Polarski and A. A. Starobinsky, Phys. Rev. D50, 6123 

(1994), astro-ph/9404061. 
[5] C. Gordon, D. Wands, B. A. Bassett, and R. Maartens, 

Phys. Rev. D63, 023506 (2000), astro-ph/0009131. 
[6] D. H. Lyth, C. Ungarelli, and D. Wands, Phys. Rev. D67, 

023503 (2003), astro-ph/0208055. 
[7] V. Bozza, M. Gasperini, M. Giovannini, and 

G. Veneziano, Phys. Lett. B543, 14 (2002), hep- 

ph/0206131. 

[8] G. Hinshaw et al. (WMAP) (2008), 0803.0732. 
[9] C. L. Reichardt et al. (2008), 0801.1491. 
[10] M. Tegmark et al. (SDSS), Phys. Rev. D74, 123507 

(2006), astro-ph/0608632. 
[11] P. Astier et al. (The SNLS), Astron. Astrophys. 447, 31 

(2006), astro-ph/0510447. 
[12] R. Bean, J. Dunkley, and E. Pierpaoli, Phys. Rev. D74, 

063503 (2006), astro-ph/0606685. 
[13] M. Bucher, K. Moodley, and N. Turok, Phys. Rev. Lett. 

87, 191301 (2001), astro-ph/0012141. 
[14] E. Komatsu et al. (WMAP) (2008), 0803.0547. 
[15] F. Feroz and M. P. Hobson (2007), 0704.3704. 
[16] F. Feroz, M. P. Hobson, and M. Bridges (2008), 

0809.3437. 

[17] M. Beltran, J. Garcia-Bellido, J. Lesgourgues, and 
A. Riazuelo, Phys. Rev. D70, 103530 (2004), astro- 
ph/0409326. 

[18] R. Keskitalo, H. Kurki-Suonio, V. Muhonen, and 



J. Valiviita, JCAP 0709, 008 (2007), astro-ph/0611917. 

[19] H. Jeffreys, Theory of Probability (Oxford University 
Press, 1961), 3rd ed. 

[20] R. Trotta, Contemp. Phys. 49, 71 (2008), 0803.4089. 

[21] M. Beltran, J. Garcia-Bellido, J. Lesgourgues, A. R. Lid- 
die, and A. Slosar, Phys. Rev. D71, 063532 (2005), astro- 
ph/0501477. 

[22] R. Trotta, Mon. Not. Roy. Astron. Soc. Lett. 375, L26 

(2007), astro-ph/0608116. 
[23] A. Lewis and A. Challinor, Phys. Rept. 429, 1 (2006), 

astro-ph/0601594. 
[24] A. Lewis, A. Challinor, and A. Lasenby, Astrophys. J. 

538, 473 (2000), astro-ph/9911177. 
[25] D. Langlois, Phys. Rev. D59, 123512 (1999), astro- 

ph/9906080. 

[26] A. R. Liddle and D. H. Lyth, Cosmological Inflation and 
Large-Scale Structure (Cosmological Inflation and Large- 
Scale Structure, by Andrew R. Liddle and David H. Lyth, 
pp. 414. ISBN 052166022X. Cambridge, UK: Cambridge 
University Press, April 2000., 2000). 

[27] N. Bartolo, S. Matarrese, and A. Riotto, Phys. Rev. D64, 
123504 (2001), astro-ph/0107502. 

[28] K. Moodley, M. Bucher, J. Dunkley, P. G. Ferreira, 
and C. Skordis, Phys. Rev. D70, 103520 (2004), astro- 
ph/0407304. 

[29] A. Lewis and S. Bridle, Phys. Rev. D66, 103511 (2002), 

astro-ph/0205436. 
[30] D. J. C. Mackay, Information Theory, Inference, 

and Learning Algorithms (Cambridge University Press, 

2003). 

[31] M. Bridges, F. Feroz, M. P. Hobson, and A. N. Lasenby 



11 



(2008), 0812.3541. 
[32] M. Beltran, J. Garcia-Bellido, J. Lesgourgues, and 

M. Viel, Phys. Rev. D72, 103515 (2005), astro- 

ph/0509209. 

[33] http : //www. rssd. esa. int/ index .php?project= 

planck 

[34] J. Dunkley, M. Bucher, P. G. Ferreira, K. Moodley, and 
C. Skordis, Mon. Not. Roy. Astron. Soc. 356, 925 (2005), 
astro-ph/0405462. 

[35] For general correlations of modes, multiple calls to 
CAMB are required. The separate CMB spectra are gen- 
erated along with the appropriate lensing-potential power 
spectrum, but lensing is not applied at this stage. The 
total unlensed CMB and lensing-potential power spectra 
are then formed and, finally, the latter is used to lens the 
former. 

[36] Examples include 1 — zaa and J^"° Pr(r) dr, where the 
latter has a uniform prior distribution. 

[37] Our individual MCMC chains did, however, fail the spec- 
tral convergence test developed in [34] , although the con- 
catenated chain did pass this test. 

[38] The value of R for which the chains may be considered 
converged is subjective; R < 0.01 is often considered to 
be well converged but reaching such a value with MCMC 
would have been computationally prohibitive in this ap- 
plication. 



